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We present measurements of differential cross sections and the analyzing powers A y , iTn, T20, 
T21, and T22 at E c . m . = 431.3 keV. In addition, an excitation function of iTn(8 c . m . = 87.8°) for 
~ (~~{ • 431.3 < Ec.m. < 2000 keV is presented. These data are compared to calculations employing realistic 

f*^ ' nucleon-nucleon interactions, both with and without three-nucleon forces. Excellent agreement with 

the tensor analyzing powers and cross section is found, while the A y and iTn data are found to be 
£S| ■ underpredicted by the calculations. 
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I. INTRODUCTION 



The three-nucleon continuum is an excellent testing ground for nuclear interactions, since both precise experiments 
and theoretical calculations are possible Q . Comparisons between measured three-nucleon observables and theoretical 
calculations provide stringent tests of the underlying nucleon-nucleon (NN) and three-nucleon (3N) interactions. 
Proton-deuteron elastic scattering at low energies offers two particular advantages: highly accurate experiments can 
t-H . be performed for a large number of observables, and theoretical calculations are simplified due to the energy being 
t-H ' below the deuteron breakup threshold and the relatively low number of partial waves required. 

Past experimental studies have focussed on the differential cross section and analyzing powers for p-d scattering 
at E c . m . = 2 MeV, just below the deuteron breakup threshold fU|. Theoretical calculations of these p-d scattering 
^ I observables using realistic NN and 3N interactions arc now routinely possible The calculations are found to 

1} ■ reproduce the differential cross section and tensor analyzing powers (TAPs) T20, T21, and T22 quite well, while the 
1— h ' vector analyzing powers (VAPs) A y and iT\\ are underpredicted by ~ 30%. The available 3N forces are found to 
have little influence, except for the J 71 " = ^ + partial wave Js],^). The situation in ncutron-deuteron scattering is 
similar ^-^|, although only the cross section and A y have been measured to date at low energies. By far the most 
r* ■ glaring difference between theory and experiment is in the VAPs, a discrepancy which has come to be known as the 
. £-i \ U A V (Q) puzzle" ||. This puzzle has been in existence for about ten years now, and the resolution is still unclear. It 
has been suggested that the discrepancy arises from inadequacies in the NN potential in the 3 Pj waves ||,|l0]]. Other 
workers have considered the role of 3N force effects |||],[ll]]. So far, no conclusive solution to the problem has been 
found, although most recent attention has been focussed on the role of the 3N interaction 12|-14|]. 



It is desirable to determine these observables at energies below E c , m . = 2 MeV, as the influence of higher partial 
waves is strongly reduced, allowing the dominant S- and P- waves to be investigated with greater clarity. Also, when 
faced with disagreement between theory and experiment, it is interesting to study the energy dependence of the 
discrepancy, since this may provide clues as to its origin. New data will place strong constraints on any proposed 
solution. There is also considerable interest in obtaining low-energy data for the determination of the p-d scattering 
lengths @. 

In this paper we present measurements of the differential cross section and analyzing powers A y , iTu, T20, T21 
and T22. f° r P"d scattering at E c . m . — 431.3 keV. These data are the first complete set of observables measured at 
such a low energy; previously only cross section Jl5|,[l6| and A y data |l(| have been obtained in this energy range. We 
also present an energy excitation function of iTu(8 c , m , = 87.8°) from E c . m , = 431.3 keV up to 2000 keV. Some of the 



data presented here have appeared previously in brief publications |lT|,[r7 18 
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II. EXPERIMENTAL PROCEDURES 



Measurements of the differential cross section and the analyzing powers A y , iTn, T20, T21, and T22 were carried out 
at the Triangle Universities Nuclear Laboratory (TUNL). The analyzing power T22 was deduced from measurements 
of T 2 o and A yy . Achieving a center-of-mass energy of 431.3 keV requires either a 650-keV proton beam incident on 
a deuteron target, or a 1300-keV deuteron beam incident on a proton target (the small corrections for energy loss in 
the target will be discussed later). The beams, targets, and detection methods are discussed below. 

A. Low-energy proton beam 

Beams of 650- keV protons were used for measurements of differential cross sections and A y . Polarized or unpolarizcd 
beams of 72-keV 1 H _ ions were produced by the TUNL atomic beam polarized ion source p9fl . These ions were then 
accelerated to 450 keV by the minitandem accelerator |2(J. The final energy of 650 keV for the now positive ions 
was achieved by the — 200-kV bias on the 107-cm-diameter scattering chamber |2l| l where the measurements were 
performed. 

Unpolarized proton beams were produced by bleeding hydrogen gas into the ionizer of the ion source. When polar- 
ized proton beams were used, a Wien filter downstream from the ion source was utilized to set the spin-quantization 
axis vertical in the laboratory. Two spin states were used, with polarizations pz ~ ±0.7. The desired hyperfine states 
of atomic hydrogen were cycled approximately every second. This technique minimizes the effects of slow changes in 
beam position, target thickness, or amplifier gain on the measured analyzing powers. Proton beam polarization was 
determined using the 6 Li(p, 3 He) 4 He reaction in a polarimeter |2^] located at the rear of the scattering chamber. The 
polarization was measured several times throughout the measurements at an incident proton energy of 450 keV by 
lowering the chamber bias voltage. The proton polarization was found to be constant within ±3% throughout the 
measurements; the systematic error in the proton polarization is estimated to be ±4%. 

The energy calibration of the proton beam produced by the minitandem and high-voltage chamber system has 
been established to ±1 keV, using the 240.0- and 340.5-keV resonances in 19 F(p, 07) and the 405.4- and 445.8-keV 
resonances in 27 Al(p, 7). The resonance energies were taken from Ref. j23j. 

B. Deuteron beam 

Polarized beams of 72-keV 2 HV ions were produced using the same atomic-beam polarized ion source [ [l9| used for 
protons. The beam was injected into the TUNL FN tandem accelerator, magnetically analyzed, and then delivered 
to a 62-cm diameter scattering chamber. 

The same Wien filter mentioned previously was used to control the spin-quantization axis on target. The spin 
axis was longitudinal for the T20 measurements, normal to the reaction plane for iT\\ and A yy ^ and 45° offset from 
longitudinal in the reaction plane for T^x- Three deuteron beam polarization states were produced with the atomic 
beam polarized ion source: a maximum positive, a maximum negative, and an unpolarized state. The TAP data were 
obtained with pzz ~ ±0.70 and pz ~ T0.25, while beams with pz ~ ±0.55 and \pzz\ < 0.05 were used for the iT\\ 
measurements. The spin states were also cycled approximately once every second. The beam polarization for the TAP 
data was determined using the 3 He(d,p) reaction in an online polarimeter located behind the scattering chamber p4[ . 
Deuteron beam vector polarization was determined online via the 12 C(d,p) reaction in a different polarimeter located 
behind the scattering chamber. The effective iTn for this reaction at Ed = 1.3 MeV has been calibrated relative to 
the 3 He(d,p) polarimeter at Ed — 12 MeV. More information on the use of the 12 C(d,p) reaction for deuteron vector 
polarimetry is available in Ref. |^5fl . For the iTn(87.8°) measurements the polarization was very stable over time, and 
was monitored periodically with the 3 He(d,p) polarimeter at Ed = 12 MeV. The absolute uncertainty in all of the 
deuteron beam polarizations is estimated to ±3%. The unpolarized deuteron and proton beams used for the 1 H(ci, d) 
cross section measurements were produced by a direct extraction negative ion source. 

The energy calibration of the beams produced by the FN tandem is determined within ±0.1% from the magnetic 
field in the 52° bending magnet, measured by an NMR magnetometer. The magnetometer was calibrated using the 
7 U(p,n) threshold (E p = 1880.443 ± 0.020 keV |§) and the 19 F(p,n) threshold {E p = 4234.3 ± 0.8 keV ||). 
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C. Targets 



The measurements were performed by bombarding thin hydrogenated or deuterated carbon foils. The targets used 
for analyzing power measurements consisted of approximately 1 x 10 18 and 1.5 x 10 18 hydrogen isotope and carbon 
atoms/cm 2 , respectively. The targets used for the relative and absolute cross section measurements consisted of 
approximately 5 x 10 17 and 8 x 10 17 hydrogen isotope and carbon atoms/cm 2 , respectively. In order to calculate 
the energy loss in the targets it is necessary to know their carbon and hydrogen isotope thickness and the stopping 
power. The carbon content was determined by comparison with carbon foils of known thickness while the hydrogen 
thickness was determined cither from the d-p elastic scattering cross section, or from the p-p elastic scattering cross 
section. For the analyzing power measurements the beam-energy loss in the targets is « 10 keV, while for the cross 
section measurements, the beam-energy loss was » 5 keV. In all cases the incident energy was adjusted so that the 
mean energy in the center of the target corresponded to E c . m . = 431.3 ± 0.8 keV. Here the quoted error includes 
the uncertainties in the incident energy and also the energy loss in the target. It should be noted that the use of 
thin targets is very important at low energies for minimizing energy loss and straggling effects. The differential cross 
section is particularly sensitive to the energy, a point which will be discussed quantitatively later in this paper. 

D. Detection methods 

The reaction products were detected using 100-/im-thick Si surface barrier detectors. Left-right symmetric detector 
configurations were utilized for all of the measurements. Count rates in the detectors were controlled by varying 
the solid angle (0.2-5 msr) and beam current (10-150 nA). In some cases, thin mylar foils (2-5 /xm) were placed 
over the detectors to either (1) increase the energy separation between proton and deuteron peaks or (2) stop heavy 
recoils resulting from 12 C(<i, 12 C)d, 12 C(d, 13 C)p, or 12 C(d, 13 N)n reactions. Measurements at far forward angles in the 
laboratory are limited by the high flux of elastic-scattered particles from carbon. For the proton-beam measurements 
in the high-voltage scattering chamber, detector signals were sent to ground potential via fiber-optic connections. 

For the cross section, T20, I21, and A yy measurements the elastic scattering yields were determined from the energy 
spectra. Sample energy spectra are shown in Figs. It should be noted that the low-energy heavy recoils seen in 
Figs. I and I can be eliminated if necessary by using the aforementioned mylar foils. The primary background under 
the p-d elastic scattering peaks at all angles is the low-energy tail from carbon elastic scattering. Pulsers were inserted 
into the detector electronics to facilitate dead-time corrections, which were typically less than 5% but always less than 
10%. Peak yields were extracted using linear least-squares fits to the background on either side of the peak. For the 
analyzing power measurements, the background subtraction and dead-time corrections were performed separately for 
each spin state. 

A different technique was used for the VAPs A y and iTn since the magnitudes of these observables are roughly 
a factor of 10 smaller than for the TAPs. The scattered deuterons and protons were detected in coincidence using 
two pairs of silicon surface barrier detectors placed at symmetric angles on either side of the incident beam. The 
angles of the detectors were set to observe either protons or deuterons in the more forward detectors in coincidence 
with deuterons or protons detected in the more backward detector on the opposite side of the beam. A major benefit 
of this method is that data can be acquired at a much higher rate because the events resulting from carbon elastic 
scattering are essentially completely eliminated. Note that a factor of 10 reduction in magnitude of the analyzing 
power requires a factor of 100 increase in the number of counts if the same relative statistical accuracy is desired. An 
additional advantage of the technique is that the background under the peak of interest is much reduced, lowering the 
possibility of false asymmetries resulting from the asymmetry in the background. Histograms of the time difference 
between the fast timing signals from each coincident pair of detectors were stored for each spin state. The time 
resolution for the coincident proton-deuteron peaks was « 8 ns, with backgrounds < 3%. A sample time spectrum is 
shown in Fig. [|. Dead-time corrections (< 3%) were determined by sending test pulses to the detector preamplifiers 
with time delays adjusted to give distinct peaks in the time spectra. Peak yields for each spin state were extracted 
using linear least-squares fits to the background on either side of the peak, and were corrected for dead time. The 
possible influence of channel (time) dependent dead-time corrections was calculated and found to be negligible. 

E. Measurements 

The incident beam energy was adjusted to yield E crn , = 431.3 keV in center of the target, except for the iTn(87.8°) 
measurements where the incident deuteron energy varied over 1.3 < Ed < 6.0 MeV. The beams were collimated to 
produce a 2 mm (horizontal) x 4 mm (vertical) beamspot on the center of the target. The 0° position for each 
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detector was determined by measuring the elastic scattering from carbon on either side of the beam for angles near 
0°. The systematic uncertainty in the angular positioning is estimated to be ±0.1°. The number of incident particles 
was determined by beam current integration. For the analyzing power data and relative cross section data the targets 
were replaced approximately every 12 hours, or when the hydrogen isotope content decreased by ps 30%. 

Relative cross sections were determined with a proton beam in the high-voltage scattering chamber. The measure- 
ments were performed using a movable pair of left-right symmetric detectors, and a fixed pair located out of the plane 
of the other detectors at 9i a b — 37.5°. The yield of protons and deuterons from p-d elastic scattering in the fixed pah- 
was used to normalize the yields in the movable pair. 

Absolute p-d cross sections were determined relative to the proton-proton elastic scattering cross section at a 
nominal proton beam energy of 2600 keV. Detector pairs were placed at 8i a b = 25° and 35° where both cross sections 
vary smoothly with angle. These angles give d-p elastic scattering measurements at c m = 82.6°, 110°, and 130°. A 
direct extraction negative ion source was used to inject beam into the tandem accelerator for these measurements. The 
beam could be switched between 1300-keV deuterons to 2600-keV protons in about one minute by only adjusting two 
parameters: (1) the magnetic field in the bending magnet between the source and the accelerator, and (2) changing 
the accelerator terminal voltage. Note that the magnetic field in the 52° analyzing magnet, which defines the beam 
energy, is unchanged throughout this procedure. The deuteron and proton beams were switched several times for 
each target, using an integrated charge of ~ 1 [iC per measurement. For these measurements accurate relative beam 
current integration is essential. Current was integrated from the target rod (biased to +100 V) and a plate at the 
rear of the chamber (suppressed by a shroud at —100 V). As a test of the technique, we also measured using the same 
procedure 197 Au(d, d) and 197 Au(p,p) elastic scattering at 6>i a b = 140°, where the cross sections are described to a 
very high accuracy by the Rutherford formula. A target consisting of 2 x 10 17 Au atoms/cm 2 was utilized for these 
measurements. The results agreed within 0.8% with the ratio expected for Rutherford scattering (« 1% corrections 
due to electron screening [Q were taken into account). 

III. DATA ANALYSIS AND RESULTS 
A. Cross section measurements 

The relative cross section data as a function of center-of-mass angle are shown in Fig. M. These data have been 
normalized to the absolute cross section measurements. 

Absolute p-d cross sections were determined at 9^ = 25° and 35° relative to the p-p scattering measurements at 
a nominal proton beam energy of 2600 keV. The p-p cross section for these energies and angles appears to be well 
understood. Our procedure is to utilize the Nijmegen energy-dependent partial- wave analysis [t29|]30[ to generate cross 
sections for our energy and angles. A thorough discussion of the data upon which this analysis is based is given in 
Ref. pi|| . The p-p cross section data for E v < 10 MeV which are included in the Nijmegen analysis generally have 
absolute uncertainties < 0.5% and are reproduced by the fit within this error. In addition there are the p-p cross 
section data of Knecht, Dahl, and Messelt (32) which include measurements at E p — 2425 keV, very close to the 
energy of the present experiment. The quoted experimental uncertainties vary in the range 0.1-0.3%. Although not 
used in the Nijmegen analysis, the Nijmegen energy-dependent fit agrees with the results of Ref. Q at 2425 keV 
within 0.4% for the angles of our interest. The systematic uncertainty in the absolute p-d cross sections is estimated 
to be ±1.1%; the contributions to this uncertainty are summarized in Table Q. The assumed p-p cross sections and 
the p-d cross sections derived from them are given in Table H 

B. Analyzing power measurements 

Analyzing powers were determined from the particle yields measured for each spin state. The yields for each spin 
state were corrected for background and dead time, and were normalized by the number of incident particles. Tests 
for false asymmetries using tensor-polarized deuteron beams were carried out by measuring 197 Au(d,d) scattering at 
@iab — 40° under identical conditions as the p-d TAP measurements. For Ed — 1.3 MeV, all of analyzing powers 
for 197 Au(d,d) are expected to be < 10~ 4 The results were consistent with zero at the level of 5 x 10~ 4 , the 

statistical uncertainty of the measurement. 

For each pair of left-right symmetric detectors, the analyzing powers were determined as follows. We first define 
the following ratios 
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where p^ and p^z are ^ ne vector and tensor polarizations for spin state i. It should be noted that Eqs. (H) and (||) are 
approximate; for the analyzing powers and polarizations encountered in this experiment the corrections are negligible. 
Values for the TAP T22 were derived from A yy and T20 using T22 = —A yy /\/?> — T2o/V&- The analyzing power data 
are shown in Figs. ^| and 0. The measured iTn(87.8°) excitation function is shown in Fig. ||. 
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IV. DISCUSSION 



These experimental results are compared to calculations utilizing the Pair-Correlated Hyperspherical Harmonic 
basis p4[ to construct the scattering wave function, and the Kohn variational principle to determine the scattering 
matrix elements [Q. The calculations were performed using the AV18 NN potential |3j| and with AV18 plus the 
3N interaction of Urbana (UR-IX) |3(|, and are shown in Figs. [| |7|, and ^|. All of the cross sections in Fig. [| have 
been divided by the cross section calculated using the AV18+UR-IX potential, so that the small differences are more 
apparent. The cross section measurements are seen to be in excellent agreement with the calculations using the 
AV18+UR-IX potential. The calculation using the AV18 potential alone differs from the measurements by 3-4%, a 
small but statistically-significant amount. As shown in Fig. ||, the TAP data are in reasonably good agreement with 
calculations using either AV18 or AV18+UR-IX. It is, however, significant to note that the AV18+UR-IX potential 
gives slightly better agreement for all of the TAPs. 

The VAP results are shown Fig. [?]. As we have previously reported the calculations using both the AV18 and 
AV18+UR-IX potentials underpredict the data by w 40%. This discrepancy has now been observed for a wide range 
of energies for both n — d and p — d scattering. The origin of this 11 A y puzzle" is not clear. Specific efforts to study 
this problem have been undertaken recently by employing new forms of 3N forces HQ. Moreover in Ref. |37|1 a new 
NN interaction constructed from chiral perturbation theory has been shown to give a better description of the VAPs 
at low energy. 

The measured energy dependence of iTn(87.8°), shown in Fig. |^, will have to be explained by any proposed solution 
to the A y puzzle. It is seen that the AV18+UR-IX potential consistently underpredicts iTu for all energies below 
the deuteron breakup threshold. Although not shown, the AV18 potential calculation lies « 10% lower than for 
AV18+UR-IX at these energies as well. We also note the the present measurements at E c . m , — 1.67 and 2.00 MeV 
are in excellent agreement with previous data . 

The comparison between the present data and theoretical calculations at E c , m , — 431.3 keV can be made more 
quantitative through the use of the x 2 parameter. For the N = 94 data points the AV18 calculation yields x 2 = 844, 
while the AV18+UR-IX calculation yields x 2 = 268. The result that x 2 /N is much greater than one even for 
the AV18+UR-IX calculation is expected since there are large discrepancies for the VAPs as well as smaller but 
statistically-significant differences in the other observables. A more detailed discussion of x 2 comparisons and also 
phase-shift analysis by x 2 minimization will be the subject of another publication. Here we comment on two systematic 
uncertainties in the experimental data and how they affect the calculated x 2 - Considerable attention was given to 
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the determination of the energy of the experiment, E c , m , = 431.3 ± 0.8 keV. In Fig. |h] we show the sizable effect on 
X 2 of varying the energy of the theoretical calculations with the AV18+UR-IX interaction. It is clearly important 
for the experimental energy to be determined as precisely as possible. We have also investigated the importance of 
the absolute normalization of the experimental cross sections. In Fig. [n] we show the effect on \ 2 °f varying the 
normalization of the cross section. Here \ 2 is calculated for the 22 cross section points. Again very significant effects 
are observed. 



V. CONCLUSIONS 



We have measured differential cross sections and the analyzing powers Ay, iT\\, T20, T21, and T22 at E CM . = 
431.3 keV, as well as an excitation function of iT u (6 c . m , = 87.8°) for 431.3 < E c . m . < 2000 keV. The 431.3-keV 
data comprise the lowest-energy complete set of cross section and analyzing power observables for nucleon-deuteron 
scattering. The cross section and TAPs are found to be in good agreement with calculations using the AV18 NN 
interaction and the Urbana-IX 3N interaction. The VAPs A y and iT\i are found to be underpredicted by ~ 40%. At 
this time it is not clear whether these discrepancies result from inadequacies in the assumed NN interaction, the 3N 
interaction, or some other source. In this context, the very low energy data presented here represent a very stringent 
constraint for any new theoretical model of the nuclear interaction. 
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FIG. 1. Charged-particle spectrum obtained at 6i a b = 30° 
with 650-keV protons incident on a carbon-deuterium target. 
The small peak near channel 380 results from H(p,p). 



TABLE I. Systematic errors in the absolute cross section 
determination. 



source 


error 


p-p cross section 


0.5% 


target uniformity 


0.5% 


beam-current integration 


0.8% 


total 


1.1% 
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FIG. 2. Charged-particle spectrum obtained at 9i a b = 25° 
with 1300-keV deuterons incident on a carbon-hydrogen tar- 
get. The peak near channel 70 also results from 1 H(ci, d), but 
the lower-energy deuteron peaks were not analyzed due to the 
larger backgrounds. The small peak near channel 270 results 
from 12 C(d,p) 13 C(3.09 MeV). The broad structures below 
channel 250 result from heavy-ion recoils from the target. 
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FIG. 3. Charged-particle spectrum obtained at 6i a b = 35° 
with 1300-keV deuterons incident on a carbon-hydrogen 
target. The small peak near channel 310 results from 
12 C(d,p) 13 C(3.09 MeV). The broad structures below chan- 
nel 250 results from heavy-ion recoils from the target. 
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FIG. 4. Proton-deuteron time-of-fiight spectrum obtained 
for 6 cm. — 87.8° with 1300-keV deuterons incident on a car- 
bon-hydrogen target. The horizontal axis calibration is ap- 
proximately 0.4 ns/channel. 
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FIG. 5. Relative proton-deuteron differential cross section 
data obtained E c . m . = 431.3 keV. The squares and circles re- 
sult from detected protons and deuterons, respectively. These 
data have been normalized to the absolute cross section mea- 
surements. 



TABLE II. Laboratory proton-deuteron cross sections with 
statistical errors determined at E d — 1294 keV together with 
the proton-proton elastic scattering cross sections assumed for 
E p = 2589 keV. 







P(P,P) 
(mb/sr) 


p(d,d) 
(mb/sr) 


P(d,p) 
(mb/sr) 


cr(6lab 


= 25.0°) 


518.9 


2240 ± 20 


926 ± 10 


a(6iab 


= 35.0°) 


493.4 




729 ±8 
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FIG. 6. The measured TAP data (circles) and theoretical 
calculations using the AV18 (dashed line) and AV18+UR-IX 
(solid line) potentials. 
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FIG. 7. The measured VAP data (circles) and theoretical 
calculations using the AV18 (dashed line) and AV18+UR-IX 
(solid line) potentials. 
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FIG. 8. The measured iTii (circles) for 9 c . m . = 87.8° plot- 
ted versus center-of-mass energy. The solid curve is a theo- 
retical calculation using the AV18+UR-IX potential. 
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FIG. 9. Ratios of cross sections to the theoretical cross 
sections calculated using the AV18+UR-IX potential. The 
experimental cross sections are shown as circles and squares 
for the relative and absolute measurements, respectively. The 
ratio of the AV18 cross section to the AV18+UR-IX cross 
section is shown by the solid line. 
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FIG. 10. Plot of x 2 versus center-of-mass energy for the 
comparison of the experimental data to theoretical calcula- 
tions at different energies. The arrow indicates the actual 
energy of the experimental data. 
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FIG. 11. Plot of x 2 versus the cross section normalization 
factor. 
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